Exact time-reversal focusing of acoustic and quantum excitations in open cavities: 

The perfect inverse filter 
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The time-reversal mirror (TRM) prescribes the reverse playback of a signal to focalize an acoustic 
excitation as a Loschmidt echo. In the quantum domain, the perfect inverse filter (PIF) processes 
this signal to ensure an exact reversion provided that the excitation originated outside the cavity 
delimited by the transducers. We show that PIF takes a simple form when the initial excitation 
is created inside this cavity. This also applies to the acoustical case, where it corrects the TRM 
and improves the design of an acoustic bazooka. We solve an open chaotic cavity modeling a 
quantum bazooka and a simple model for a Helmholtz resonator, showing that the PIF becomes 
decisive to compensate the group velocities involved in a highly localized excitation and to achieve 
subwavelength resolution. 

PACS numbers: 03.65.Wj; 03.67.-a; 43.60.-c 
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INTRODUCTION 

The time-reversal mirror (TRM) technique was devel- 
oped by Mathias Fink and collab. |1] in order to refocus 
an ultrasonic excitation diffusing away through a hetero- 
geneous medium. The remarkable robustness of this tech- 
nique against inhomogeneities and imperfections in the 
procedure enabled novel applications in medical physics 
[2(], oceanography [3( and telecommunications A]. The 
acoustic signals emerging from a localized source point 
are recorded by an array of transducers. Ideally, this 
array splits the system in two portions: the internal cav- 
ity, where the focusing takes place, and the outer re- 
gion, towards where the excitation finally spreads away. 
The transducers then play-back each registry in the in- 
verted time sequence. Surprisingly, this simple prescrip- 
tion allows recovering the initial excitation with high pre- 
cision. This constitutes an alternative form of achieving 
a Loschmidt Echo [5|, i.e. the time reversal of an ex- 
citation distributed in a complex medium. To ensure 
the exact reversed dynamics, we exploited the close cor- 
respondence between the Schrodinger and the classical 
wave equations @, 0] by introducing a formal prescrip- 
tion denoted as the perfect inverse filter (PIF) [8]. This 
yields the precise injection function that produces the 
targeted dynamics inside the cavity. In particular, the 
PIF prescribes how to process the recorded signal in or- 
der to achieve the exact refocusing of whatever excitation 
exists in the cavity, provided that it was originated as a 
wave packet incoming from the outer region. A poten- 
tial application is the quantum bazooka, analogous to the 
acoustic bazooka conceived by M. Fink and collaborators 
0. When a waveguide injects a well-defined wave packet 
into a chaotic cavity, it is reflected as a persistent noisy 
signal. The time-reversed counterpart allows to store in 
the cavity the weak signal injected through a long period 
of time until it emerges as a burst of excitation. 

To further optimize the spatial contrast of a localized 



acoustic excitation in cases where TRM results imperfect, 
Fink and collab. proposed the spatio-temporal inverse fil- 
ter (STIF) [10]. It requires the inversion of a matrix built 
from the Fourier transforms of the direct wavefield prop- 
agators. These give the response at the transducers to a 
Dirac delta signal in each control point within the focal 
region. Unlike the TRM, the STIF results in a highly 
invasive technique. Minimally invasive STIF (miSTIF), 
a recent proposal by Vignon et al. [Ill ], achieves the 
same focusing quality as STIF without resorting to con- 
trol points. All the information in the wavefield propa- 
gator is deduced from the backscattering signals of the 
transducers and the technique results no more invasive 
than TRM. Both the miSTIF and the PIF involve the 
inversion of a propagator matrix in the reduced basis of 
the transducers. However, this may result in numerical 
instabilities that require a delicate regularization proce- 
dure to deal with the Green's function divergencies. 

In this work, we assume a natural separation between 
the incoming and outgoing excitation components, as 
when the source is internal, to develop a PIF procedure 
based on the detection of the escaping wave. This in- 
ternal PIF results very simple and intrinsically stable. 
Moreover, it reveals a property not obvious before: the 
filter that corrects the TRM does not depend on the in- 
ternal scattering, but only on the group velocities of the 
various propagating modes in the far field region. This is 
general enough to cover, either exactly or as an approx- 
imation, a wide range of TRM set-ups that would bene- 
fit from a PIF. We may mention subwavelength focusing 
[IH , lithotripsy in dispersive media , the quantification 
of nonlinear media jl 31 ] , coherent control of optical excita- 
tion energies in nanoplasmonic systems 14[ and acoustic 
bazookas @. To this list we may add the development 
of novel coherent control strategies that seek the local 
injection of designed quantum excitations. As discussed 
below, this involves fields as diverse as nanoelectrome- 
chanics, Bose-Einstein condensates and NMR. 
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We start by considering a quantum system where the 
corresponding propagators are well known (IBj . Then, 
we will show that the resulting PIF prescription remains 
valid to describe sound propagation. A numerical con- 
trast with the TRM procedure in two situations of phys- 
ical interest, the quantum bazooka configuration and a 
simplified model for a Helmholtz resonator coupled to 
a grated waveguide, confirms the relevance of the PIF 
correction. 



TIME-REVERSAL IN A CAVITY 

We consider the situation in which, at the initial time, 
the excitation is concentrated inside the cavity. From 
to = until a registration time t R , the discrete wave 
function 



E, 



mG?„(t)!Mo), o < t < t R , (i) 



is detected and recorded as it escapes through the fron- 
tier. In the present letter we consider a cavity, cither intc- 
grable or chaotic, open through a Id channel. Therefore, 
its frontier is controlled by a single transducer at x s = sa, 
with s integer and a the distance between neighbors sites 
in a lattice. Although apparently a restriction, this sim- 
ple model faithfully represents the essential underlying 
physics. A straightforward generalization requires a ma- 
trix formulation and will be presented elsewhere. 

Here, n (t) is the retarded Green's function relating 
the amplitude response at x m when a Dirac delta signal 
8(t) is injected at x n inside the cavity. It satisfies the 
inhomogencous Schrodinger equation 

M§- t G*, n (t) - H m ,iG* n (t) = S m , n S(t), (2) 

with H rn i the Hamiltonian matrix elements. Assuming 
that there are no localized states, the registration time 
t R is taken long enough in order to ensure that all exci- 
tations have left the cavity. 

As we have shown in ref.Q, the time reversion for a 
wave packet arriving at the cavity from the outer region 
is performed, in the energy domain, through the injection 



Xs{e) 
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Notice that in the above expression, 



(3) 



(4) 



accounts for the time reversed complete evolution of the 
wave packet, i.e. both the incoming and outgoing com- 
ponents are required in order to compute ip r s ev (t). In the 
present work we want to reverse the signal that is pro- 
duced inside the cavity. Then, we deal with the building 



up of the complete evolution from the knowledge of ip s (t) 
at times < t < t R . On the assumption that the injec- 
tion function is known, the reversed propagation at the 
transducer would be 

ipr(t) = r s (2t R -t) (5) 

= ih f Gf s (t-t') x (x S7 t')dt\ (6) 
Jo 

for times t R <t < 2t R before the focalization. Since the 
PIF prescription achieves exact reversion inside the cav- 
ity, the evolution at subsequent times (t > 2t R ) can be 
conceived as the wave packet that starts at the focaliza- 
tion time with ip n (2tR) = ipn(0). Thus, for subsequent 
times we expect 



^ s CV (t)=J2^G?,n(t-2tRWn(0). 



(7) 



A comparison between eqs.Q} and ([7]) shows that for 
the most simple case in which all coefficients ipn(0) are 
real, both evolutions become identical. The same evo- 
lution but sign changed should be obtained when they 
are purely imaginary. In a more general situation where 
the coefficients present different phases we start with the 
Fourier transform 
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\p* s (2t R -t)e iet/h dt+ 



(8) 



+ ihJ2 G« n (t-2t R W n (0)e^ h dt (9) 

„ J2tn 



Je2t R /h 



OOC(o) 



(10) 



where the second term within the brackets can be inter- 
preted as the unknown evolution for subsequent times. 
The e argument and the superscript in G^ n (e) are omit- 
ted for a compact notation henceforth. At this point, we 
account for the hopping terms V s , s -i and V^-i^ (with 
x s -i inside de cavity) connecting the cavity with the 
outer region. Hence, we use the Dyson equation [l6| to 
rewrite 



G s 



G?l + G s 



s fs.s-l^s-i „, 



(11) 



where G^ follows eq.((2]) for the unperturbed Hamilto- 
nian with V s ^ s -i — V s _i )S — 0. Since s and n are at dis- 
connected subsystems G s °n = 0. Besides, since the cavity 
subsystem is closed and non absorbing, V s , s ~iG < f} 1 n is 
a real number. Therefore, we obtain for the unknown 
evolution 

mJ2G s , n r n (0) = -^r s (e). (12a) 

Finally, we nse this last expression to rewrite eq. ffTO)) as 

G s 



Vs ev (e) 



Je2t R /h 



1 



We). 



(13) 
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Since we have the complete 4> r s ev {£), the injection function 
according to the PIF formula is given by 



Me) 
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r 



Je2t 



^T(e) (14) 
«/*Im(G-i)#(e). (15) 



This is the PIF prescription for the case where the initial 
state is an excitation inside the cavity. As we shall see, 
the imaginary component of the inverse in G StS is closely 
related with the group velocity of the scattered waves 
in the outer region. Therefore, contrary to the external 
source condition presented in ref. Q , the correction im- 
posed by the PIF in internal source condition does not 
depend on the structure of the cavity. 

In this Letter, we used the one-dimensional picture 
for simplicity. Following the same reasoning as before, 
a multichannel generalization for the PIF procedure im- 
plies minor details and will be described elsewhere. In 
such case, one should be able to compute the relevant 
Green's function components linking the transducers in 
the boundary through a matrix continued fraction algo- 
rithm 17]. An alternative, more general formulation for 
the PIF may be obtained based on the fundamental opti- 
cal theorem 18j. In words, it states that the continuous 



density of states at a given site xq is built upon the es- 
cape rates through the boundaries of the system. If we 
choose a semi-infinite chain and there is no escape to the 
left of site xq, the excitation can only escape to the right 
of the frontier x s with a rate ImS s /?i. In equations: 



ImG 0l o = G , s lmE s G 



s,o- 



(16) 



Multiplying both sides by a pulse £(e) and the ift factor 
we find 

ifi (G ,o - G$ )0 ) m = 2iG , s r s (-i7i)G s * £(e), (17) 

with r s = — Im£ s . Since the outgoing wave originated 
at the focal point writes 



(18) 



we re-obtain the PIF formula by identifying the second 
term in the left hand of eq. (|17[) as the incoming wave 
that produced the pulse in the focal point. The complete 
evolution, identical to the reversed one, results 



(19) 
(20) 



The term between the square brackets is the excitation 
that must be injected at x s to reproduce the original sig- 
nal in the cavity and coincides with the PIF prescription 
obtained in eg. (ITS]). 



PIF in a quantum bazooka configuration 



In order to evaluate the eq. (|T5t we propose a model for 
a quantum bazooka device composed by a stadium bil- 
liard coupled with a one-dimensional waveguide as shown 
in the top of fig. ([1]). The discretization in the Schrodinger 
equation gives the tight-binding structure of the Hamil- 
tonian 



(21) 

where Hb is the Hamiltonian of the billard in a discrete 
basis and Hbc the coupling with the waveguide. In the 
waveguide, Ej denotes the energy at site Xj — ja and 
Vjj+i is the hopping amplitude between sites Xj and 
Xj+\. Supposing a single semi-infinite waveguide, the 
cavity results delimited by a single transducer at site x s . 
The Green's function is obtained through the continued 
fraction technique as 



1 



£-£ s -£ in (e)-£ out (e)' 



(22) 



where S; n (e) and S out (e) are the self-energy corrections 
due to the presence of sites inside and outside the cavity 
respectively [3] . Notice that the presence of the billard 
is included through Ei„(e). 

Decimation on the cavity gives the Si n (e) contribution 
as a continued fraction composed by the hoppings and the 
site energies. In absence of magnetic fields or dissipation 
those parameters are all real numbers. Therefore, 



Sin(e) = 



V a , a -iV a -i, 



E s -i 



s-2K,_2,s-l 



(23) 



e - E s -2 - ■ ■ 



Vi.oVb, 
e-E -. 



is also a real function regardless the details of H B . 

On the other hand, in the limit where the number of 
sites increases indefinitely, the homogeneous outer region, 
a linear chain with site energies E and hoppings V, con- 
tributes to the self energy as a complex function [l9j : 



t (e)=A( £ )-ir( £ ). 



(24) 



The relationship between T(e) and the group velocity can 
be found through the dispersion relation for the asymp- 
totic waves: 



£fc = E Q — 2V cos(ka), 

__ 1 de k _ IV a 
Vg ~ h~dk ~ ~JT 



(25) 



sin(fca) = i> m ax sin(fca) (26) 




We choose E Q — 2V so that the differential form of 
the Schrodinger equation for a particle with mass m is 
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FIG. 1: (color online). Focalization at the focal point for 
TRM (red dashed) and PIF (blue dotted) procedures con- 
trasted with the exact reversal (solid). PIF and exact curves 
coincide for times close to 2tn. Top: basic scheme of a quan- 
tum bazooka device with a localized excitation inside. Left 
inset: spatial contrast in the recovered signal. 



obtained as the limit h 2 /2Va 2 — > m when V — >• 00 
and a 2 — > 0. Since only the imaginary component from 
^out(^) is required, eq. (jT5]l writes 

x ™(e) = |e ie2t «/ R r(e)V;(e). (28) 

The complex exponential serves to define the origin of 
time, and hence we can neglect it. Using the definition 
of r(e) given above we get 



= 11- 



2V 



2V 



WW 



-X 



TRM 



(e). 



(29) 
(30) 



The exact reversal requires the injection of a TRM sig- 
nal filtered by the group velocity of the scattered waves. 
When the initial state is a local excitation composed by 
a few sites inside the cavity, the Fourier transform of the 
detected signal ip s (e) will cover the whole energy band. 
The factor in the PIF procedure shows that in disper- 
sive media the correction is effective near the band edges 
where the group velocity becomes negligible. 

In the numerical solution, the propagation of an initial 
gaussian wave-packet with enegy (k a) 2 V is recorded by 
a single transducer placed in the waveguide consisting of 
a single propagating mode (in practice, this waveguide is 
well represented by a one-dimensional chain). The time 
evolution was performed through a second order Trotter- 
Suzuki algorithm Fig.((TJ) shows the recovered signal 
at the focal point, i.e. the centroid of the original wave- 
packet. The focalization functions for both methods are 



close to the exact reversed propagation. However, a slight 
broadening of the focalized wave-packet is observed in 
the TRM. In this case, the differences between PIF and 
TRM procedures are small because the initial excitation 
is mainly composed by states whose group velocity in the 
outgoing channel remains almost constant. 

PIF in a classical Helmholtz resonator 

The PIF procedure for a classical wave equation can 
be readily deduced using our Green's function strategy 
by considering its finite difference version. This is best 
visualized as a system of coupled oscillators [FoJ. By in- 
cluding an inhomogeneity, it becomes in a simple model 
for a Helmholtz resonator coupled to an acoustic waveg- 
uide. Here, the lowest frequency ujq of the A/2 mode in 
the resonator is represented by a single mass too and its 
corresponding spring. The waveguide is modelled by a 
semi-infinite chain of identical masses to = am,Q placed 
at the equilibrium points x n — na (with a the lattice 
constant and n a positive integer). The nearest neighbor 
spring constant is K = muj 2 . The equations of motion 
for the corresponding displacements u n can be written in 
a matrix form which, in the frequency domain, reads: 



/ 



V 



~,2 



J 2 - 2uj 2 



u(w) = 0, 



where W 2 . 



(31) 

uj 2 + auj 2 . Here, represents a 

momentum-displacement response function, which is the 
resolvent of the dynamical matrix. Any diagonal element 
has a simple expression in terms of continued fractions. 
In particular, a far-field transducer is placed at the site 
x s on the waveguide, where the response function is 



D s . s {uj) 



1 



- 2lu 2 - A l (lu) - [A r (lu) - iw77M] ' 



(32) 

Here, the mean-field frequency 2u> 2 appears shifted by the 
dynamical effect from the oscillators at both sides of x s . 
The imaginary shift indicates that excitation components 
of different frequencies would eventually escape through 
the waveguide at the right with group velocities 



(33) 



The propagation of an excitation originated in the res- 
onator xo is detected in x s inside the waveguide. The 
recorded signal presents a strong component in ujq, which 
is the "carrier" frequency as can be appreciated from the 
density of states at the resonator's site. The displacement 
at the transducer, resulting from an initial excitation £0 
in the resonator, can be expressed as 



u.(t) = G s ,o(i)£o, t > 0, 



(34) 
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where G St o(t) is the Green's function describing the 

displacement-displacement response. In general, the 

connection between the Green's function and the 
momentum-displacement response is given by 



(35) 



As before, the complete evolution u s (t) in the transducer 
can be conceived as the forward evolution corresponding 
to the positive times and a backward evolution account- 
ing for the negative ones. In this sense, 



u s (t) 



u a (-t), t < 0, 
u s (t), t>0, 



(36) 



and the injection that produces the desired reversion can 
be obtained as 



5u a (u)) 



G s . s (uj) ' 



(37) 



According to cq. (j3~4")) . the complete evolution writes in 
the frequency domain as 

u*(uj) = u* B (u) + u B (u) (38) 
= [G* fi (oj) + G sfi (oj)]Au . (39) 

Here again, the key tool is the Dyson equation connecting 
the two subspaces delimited by the transducer. We can 
write 



G Sj0 (w) = 



(40) 



and the PIF injection rewrites in terms of the partial 
evolution corresponding to the detected signal 



1 



1 



1 



G s , s (cj) 
1 



G* 



i 



D a , a (u) 
r)(u))u*(ui). 



(41) 

(42) 
(43) 



As in the quantum case, we denote 2uj x u*(uj) as the 
Fourier transform 5uJ RM (w) of the injection function in 
the TRM protocol. In consequence, the perfect time re- 
versal is obtained only once a further filter v g ((jj)/v m ax 
is applied. Hence, the PIF formula for internal source in 
the acoustic case is 



Su 



TRM 



(44) 



Notably, the prescription remains exactly the same as 
that of the quantum version (see eq . ([30)) ) . This implies 
that the effectiveness of the filter will depend on the 
structure of the waveguide and the initial wave-packet, 
regardless the details of the cavity. However, for cases in 



which the group velocity in the free space is constant, it 
is clear that 5u™{uj) = SuJ rm (uj). 

A numerical simulation of the reconstruction of the 
initial excitation was performed using the Pair Parti- 
tioning method (PPM) 21], which yields the complete 
dynamics by alternating among the evolutions of pairs 
of coupled masses. In a similar way as the Trotter al- 
gorithm, PPM approximates the actual evolution de- 
termined by the Hamiltonian Hi + H2 during a small 
time step St as a sequence of unitary transformations 
U [(#! + H 2 )St] ~ U [Hrft] U [H 2 St] and results in a per- 
fectly reversible algorithm. The calculated dynamics is 
best depicted, as in fig. @ , by analyzing the local energy 
which avoids the fast fluctuations shown by displacement 
and momentum. Here, the left panel shows (in a log 
scale) that the recovered resonator's local energy coin- 
cides with the exact reversal over all the relevant ranges. 
Since in this model the waveguide has a cut-off frequency, 
as would be the case in a grated waveguide, the PIF filter 
improves the TRM focusing when it comes to reproduce 
the low intensity signals. This is particularly evident in 
the reproduction of the time-reversed survival collapse. 
This is a sudden dip in the local energy resulting from 
the interference between the excitation surviving in the 
resonator and that returning from the waveguide. This 
surprising phcnomcnom was originally described in the 
context of quantum spin channels [23]. The perfect con- 
trast of the time reversed signal provided by PIF is evi- 
denced in the right panel by the exact cancellation of dis- 
placements except for the resonator. It is also interesting 
to notice that both procedures produce a phantom sig- 
nal outside the "silence region" . While the TRM has an 
evident imperfection in the localization of the signal, the 
PIF procedure yields this absolute cancellation even out- 
side the cavity region defined by the far field transducer 
at x s . The resulting localization, which corresponds to 
A/ 2 of the carrier signal, was only enabled by filtering out 
the band edge components in the emitted signal. From 
this perspective, the PIF procedure contributes to the 
goal of achieving focusing beyond the diffraction limit 



CONCLUSIONS 

We presented an expression of the PIF procedure which 
results particularly simple when the excitation is gener- 
ated by an internal source. We observed that, contrary 
to what happens when the source is external, the pre- 
scription does not involve internal details of the cavity 
but only the simpler information about the propagation 
in the outer region. Hence, this filter applies to two phys- 
ically relevant situations: 1) When the excitation is actu- 
ally originated in the interior of the cavity. In this case, 
one obtains a perfect reversal of the wave function for all 
the times after the source has been turned-off. 2) When 
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FIG. 2: (color online), (a) Balls and springs model for a 
Helmholtz resonator coupled to a waveguide, (b) Focalization 
signal in the resonator. Left panel: local energy recovered 
as function of time. Both methods, the TRM (red dashed) 
and PIF (blue dotted) are contrasted with the exact reversal 
(solid). PIF and exact curves coincide for the whole temporal 
range. Right Panel: Displacements at the focalization time. 



one uses an external source in a situation that allows for 
a clear separation between incoming and outgoing waves. 
In that case, one can use the recording of the outgoing 
wave to perfectly reverse the whole excursion of the ex- 
citation through the cavity. This condition is achieved 
when the boundaries are placed far enough from the re- 
verberant region, as in the quantum and acoustic bazooka 
devices. An interesting consequence of our result in the 
acoustical case is that TRM produces a perfect time re- 
versal within the cavity, provided that the propagation 
beyond its contour is free of further reverberances. 

A numerical assesment of the reversion fidelity shows 
that PIF constitutes a notable improvement over TRM 
prescription in cases where the energy (frequency) de- 
pendence has a non-linear dispersion relation, as typical 
excitations described by the Schrodinger equation or, in 
acoustic systems, when the escape velocity is not con- 
stant. This occurs in the presence of multiple propagat- 
ing channels or when collisions outside the cavity have a 
relevant contribution as in grated waveguides. 

The PIF prescription becomes an analytical tool to de- 
sign specific excitations well beyond the discussed TRM 
context, thus providing an alternative strategy for coher- 
ent control. In particular, it allows the design of quantum 
bazookas that shoot wave packet excitations. Feasible 
implementations of this concept include the local genera- 
tion of vibrational waves in nanoelectromechanical struc- 
tures [23j . There, the excitation is injected through the 
coupling of a Cooper-pair box to a nanocantilever. In 
Bose-Einstein condensates confined to an optical lattice 
[24| . there is a well-defined set of quantum states with 
tunable couplings. Thus, the generation of macroscopic 
wavefunctions would benefit from our simple and consis- 
tent PIF prescription. Last but not least, in NMR, the 



injection of targeted spin- wave packets in chains of inter- 
acting spins |25j | is possible through the local injection of 
the polarization stored in a rare 13 C nucleus |26| . 
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